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ABSTRACT 

High resolution mid-infrared (MIR) images of the central region of NGC1068 have 
been obtained with VISIR, the multi-mode MIR instrument recently installed at the 
ESO/VLT on Paranal. A map of the emission at 12.8 urn with increased sensitivity 
over the central 8x8 arcsec 2 area is discussed. It shows a central core (unresolved along 
the E-W direction) and an extended emission which draws a spiral pattern similar to 
that observed on near-infrared images. Patches of MIR emission can be detected up to 
a distance of 4 arcsec from the core. The deconvolved 12.8 urn map is fully consistent 
with previous high-resolution MIR observations. It highlights the structure of the 
extended emission, already seen on the un-deconvolved image, and allows to identify 
a set of mid-infrared sources: 7 in the NE quadrant and 5 in the SW quadrant. The 
MIR emission map is compared with those obtained at comparable angular resolution 
in the near-infrared and in the [OIII] line emission. The very good correlation between 
the VISIR map and the HST optical map supports the idea that the MIR emission 
not associated with the torus arises from dust associated with the narrow line region 
clouds. The N-S extension of the MIR core (0.44 arcsec) is then probably simply due 
to the mixing of the MIR emission from the dusty torus and the MIR emission from 
NLR cloud B, located only 0.1 arcsec to the North. 
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1 NGC1068, A PRIVILEGED TARGET FOR 
THE VLT/VISIR SCIENCE 
DEMONSTRATION PROGRAM 

Thanks to its brightness and proximity (14.4 Mpc; 

1 arcsec=70 pc), the active galactic nucleus (AGN) in 

NGC1068 has been studied in more details than any other 

AGN. It has played a key-role in the emergence of the 

unified AGN model: Antonucci & Miller (1985) showed 

that the polarised optical spectrum of this type 2 AGN was 

identical to that of a typel. To understand this discovery, 

the authors proposed that a type 1 AGN was embedded 

within an optically thick dusty torus, hiding the central 

engine and the broad line region (BLR) from direct view 

along specific lines of sight. In parallel, the dusty torus was 

invoked to explain the type-dependent infrared (IR) bump 

observed in the spectral energy distribution of AGN and the 

conical shape of the Narrow Line Region (NLR). The work 
i 1 1 ' 1 1 1 h i 

of modellers (Krolik & Beselman 1986; Pier fe KrorikHlQa" 1 
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Granato ^TJanesej |W94 lEfstathiou fc Rowan-RobinsOx 
1995; iGranato et alJ 119971) gave support to this interpre- 
tation and a unified AGN picture involving a dusty torus 
became commonly accepted in the community. 
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Still, after several years of search, the existence of the 
torus has not been established in a decisive and direct way. 
Therefore, NGC1068 remains a natural priority target for 
the latest generation of high spatial resolution IR instru- 
ments. The high quality adaptive optics im ages obtained 
with NACO at the VLT jRouan et alJl2004h in the K, L 
and M bands (0.06 arcsec resolution in K) did not signif- 
icantly resolve the central core in the E-W direction (see 
Sect. 14. it . Recently, NGC1068 became th e first extrag alac- 
tic object to be observ ed with the VLTI iWittkowski et all 
l2004(jaffe et alJl2004) : for the first time, spatially resolved 
informations on the central emission core were obtained. 
However, the visibility measurements collected so far remain 
too scarce for a decisive interpretation and the existence of 
the torus in this AGN is still under discussion. 

Besides the core emission and the reality of the torus, 
the existence of an IR extended emissio n region around th e 
central engine is kn o wn since the works oflBock et alHl99Slh 
lAlloin et all feOOd) . iBock et all fcOOOl) ar i d iTomono eTal 
2001) i n the MIR, andlRouan et alJ Jl99sh.lMarco fc Alloin 
l200Cli. iThompson et all (l200llh iRouan et alJ (120041) and 
lGratodoiii^t^lTl2005F in the NIR. On a large scale (up to 
5 arcsec), the emiss ion is surprisingly coinc ident with the ra- 
dio jet observed bv lGallimore et al.Hl996T) . At smaller scale 



(0.5 arcsec), it appears to be coincident with the NLR cloudy 
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NGC1068 Jy/arcsec* 
VLT/VISIR, I2.8|im 11518 




Figure 1. For all the images: The X- and Y-axis a re in arcsec, North is u p and East is left, and the (0 arcsec;0 arcsec) position corresponds 
to the location of the central engine as derived in Gall iano et alj 12003). coincident with the MIR core. Images 1 to 5 were recorded by 
VISIR at the VLT. The colour scale represents either the continuum flux density surface brightness in Jy arcsec -2 (all images but 4 and 8) 
or the line flux surface brightness in 10 — 14 ergs -1 cm -2 arcsec -2 (images 4 and 8). The circles are drawn at the positions listed in tablelTl 
and have a diameter of 0.6 arcsec. (1): 12.8 [am un-deconvolved image of the central region of NGC1068 by VISIR. Contour levels are 0.07, 
0.08, 0.13, 0.30, 0.76, 1.7, 4.4, 7.2, 13.1, 23.1, 35.5, 55.1 Jy arcsec -2 . The yellow tick-marks represent the position angles measured through 
ellipse fitting at the distance from the peak corresponding to the centre of the segment. In the top-right insert, a sketch of the image has 
been drawn identifying and naming the MIR emission sources referred to in the text. (2) central part of the 12.8 urn image shown in (1). 
(3): deconvolv ed 12.8 um image. (4 ): deconvolved continuum subtra cted [Nelll image. ( 5): 18.7 urn filter image. (6) Subaru/MIRTOS 
12.3 um image iTomono et al.l200ll). (1) VLT /N ACQ 3.5 um image bv lRouan etafl 120041). (8) HST/WF PC2 F502N filter ([OIII] 5007A) 
image (cloud names after lEvans et al"lll99ll) . (9) Subaru/MIRTOS 18.5 um image iTomono et"aill200 J) . 
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structure observed by HST iMacche tto et all Il994|) . usin g 
the relative positioning discussed bv lGaliianoet"aTl 112003ft . 

The physical processes at the origin of the extended MIR 
emission still need to be clarified. Even if the fact that sili- 
cate absorption is detected all over the MIR emitting region 
llTomono et aljl200lt iGalliano et al.ll2003l) supports an in- 
terpretation in terms of warm/hot dust, the coincidence of 
the MIR emission with the radio emission, on large scale, is 
not fully understood. 

VISIR, the new MIR instrument at the VLT, received 
first light in May 2004. NGC1068 was an obvious target for 
the Science Verification Program. Preliminary results based 
on these observations are presented in this Letter. 



2 OBSERVATIONS AND DATA REDUCTION 

VISIR 1 is the recently commissioned MIR instrument in- 
stalled at the Cassegrain focus of Melipal at the VLT. It 
is composed of an imager and a long-slit spectrometer. In 
its imaging mode, a large set of filters is available, cover- 
ing both the N- and Q-bands. The possible pixel scales are 
0.075 arcsec pixel -1 and 0.127 arcsec pixel -1 respectively. 
These sizes are small enough for efficient deconvolution. 

NGC1068 was observed between 2004 November 24 
and November 30 through four filters: the 12.8 (xm fil- 
ter (A=12.81 |xm, AA=0.21 |xm), the reference filters 
for [Nell], 12.3 U-m (A=12.27 fim, AA=0.18ujn) and 
13|xm (A=13.04u.m, AA=0.22ixm), and the 18.7 |xm fil- 
ter (A=18.72 urn, AA=0.88 ixm). For the observations, the 
nodding was set parallel to the chopping, with throws of 
12.5 arcsec in the N-S direction. The total int egration time 
on so urce was 1500 s. Standard stars from ICohen et alJ 
( 1999) were observed before and after each target exposure, 
for PSF determination and flux calibration. The weather 
conditions were good, and the images are close to the 
diffraction-limit . 

Dedicated IDL routines have been written for the reduc- 
tion of the raw data. This basically consisted in shifting and 
stacking the individual frames, each frame corresponding to 
one chopping position. The precision of the flux calibration is 
limited by variations of the conversion factor due to changes 
in the transparency of the atmosphere. The final uncertainty 
on the flux of a high signal-to-noise source is of the order 
of 20%. The deconvolution was perfo rmed using the mu l- 
tiscale maximum entropy algorithm bv lStarck et alJ {2002). 
The subtraction of the continuum from the non- deconvolved 
12.8 u.m image was performed as follows: a 12.8 u.m contin- 
uum image was interpolated from the 12.3 U-m and 13 u.m im- 
ages and then subtracted from the 12.8 u.m image (all images 
matching the same resolution). For this procedure, the rel- 
ative positioning of the maps was performed according to 
the N E heart-shaped structure, pointed out bv lAlloin et alJ 

feoooD . 

3 THE DATA 

The signal-to-noise ratio of the observed PSF only al- 
lowed, among the four images obtained, the 12.8 ixm im- 

1 http:/ /www. eso.org/instruments/visir/ 



Table 1. Positions and 12.8 |xm flux densities of the sources in 
the central region of NGC1068 (from the VISIR image) 



source 


Act 


AS 


error 


Flux 




[arcsec] 


[arcsec] 


[arcsec] 


[Jy] 


nd NE7 


3.07 


2.40 


0.08 


0.05±0.04* 


nd NE6 


2.47 


2.40 


0.08 


0.05±0.01* 


nd NE5 


2.02 


2.70 


0.08 


0.07±0.05* 


nd NE4 


0.67 


1.35 


0.02 


l.±0.35 


d NE3 


0.15 


0.82 


0.02 


3.8±1.1 


d NE2 


0.47 


0.43 


0.02 


1.45±0.3 


d NE1 


-0.05 


0.49 


0.02 


11.2±2 


d Core 


0.0 


0.0 


0.0 


25±5 


d SW1 


0.14 


-0.52 


0.02 


1.2±0.2 


nd SW2 


-0.24 


-1.14 


0.02 


0.42±0.15 


nd g W3 


-1.12 


-1.94 


0.08 


0.08±0.2* 


nd g W4 


-1.51 


-3.38 


0.08 


0.04±0.03* 


nd g W5 


-0.38 


-4.05 


0.08 


0.23±0.1* 



The flux densities are measured through apertures of 0.6 arcsec di- 
ameter, centred at the position given in the table. The errors take 
into account the uncertainty in the photometric calibration. For 
the weakest more extended sources marked with a *, the very 
large quoted errors reflect the fact that the source borders are 
difficult to define, and hence a measurement through a circular 
aperture is uncertain. Still for these sources, the detection is clear, 
at a level of more than 6<r above the background noise. 
( d ): measurements performed on the deconvolved image 
(""): measurements performed on the non-deconvolved image 

age to be deconvolved. The reduced data are presented in 
Fig Q (1) 8x8 arcsec 2 central region of the 12.8 ixm non- 
deconvolved image with identification of the sources, (2) 
12.8 (i.m non- deconvolved image (same as (1), but only the 
4x4 arcsec 2 central region), (3) 12.8 |xm deconvolved image, 
(4) [Nell] continuum subtracted image, (5) 18.7 ^.m image. 

The final resolution on the images can be measured 
by the E-W extension of the core, which is known to be 
unresolved down to 0.1 arcsec (see Sect. 0. The measured 
full widths at half maximum (FWHM) are respectively 
0.29 arcsec for images (1) and (2), 0.19 arcsec for image (3), 
0.26 arcsec for image (4), 0.66 arcsec for image (5). 

In a 4 arcsec diameter aperture centred on the core, we 
measure the following fluxes: 43±8 Jy in the 12.8 ixm fil- 
ter a nd 45±9 Jy i n the 18.7 ixm filter. In the same aper- 
ture. iBock et all E oOO) found 38±4Jy at 12.5 u.m, and 
iTomono et all <200lh found 33.3±4Jy at 12.3 u.m and 
46.6±4 Jy at 18.5 ixm. The larger flux value in the 
12.8 [Lia filter can be attributed to the presence of the 
[Nell] line, which contributes to about 10% of the to- 
tal flux. From the continuum subtracted [Nell] image, we 
find a [Nell] line flux of 2.1xl0 -12 ergs -1 cm -2 through 
the same 4 arcsec diameter aperture. This flux is consis- 
tent with the ISO/SWS [Nell] 12.8 um flux: it represents 
30% of 7x 10 -12 erg s -1 cm -2 measured bv lLutz et all feOOCft 
through a much larger aperture (14x20 arcsec 2 ). 

On the 12.8 |xm image (Fig. 02), a diffraction ring 
can clearly be seen with a radius of about 1.2 arcsec. A 
residual of this ring structure remains on the continuum 
subtracted [Nell] image shown in Fig. 04. The complex 
morphology revealed by the VISIR images is consistent 
with previous observations. The image shown in Fig. 01 
details the very extended emission already observed by 
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lAlloin et ail fcOOOl) . On the 4 x4arcsec 2 region around the 
core, the deconvolved 12.8 [im VISIR image (Fig. 03) is 
strikingly simi l ar to the Subaru/MIRTOS 12.3 [im map by 
ITomono et alJ (1200 ll) . shown in Fig. 06. As we did not 
have a PSF with high enough signal-to-noise ratio to decon- 
volve the VISIR 18.7 um image, we display in stead the Sub- 
aru/M IRTOS deconvolved 18.5 urn image bv lTomono et all 
fcOOll) in Fig. 09. 

In the top-right insert of Fig. 01, we present the iden- 
tification of 13 sources, including the core. Sources SW3, 
SW4, SW5, NE4, NE5, NE6 and NE7 were identified on 
the non- deconvolved 12.8 (J.m image, at a detection level 
better than 6a above the background noise (6<r for source 
SW5). Sources SW1, NE1, NE2, NE3 were identified on 
the deconvolved image, and are clearly real since they can 
also be seen on the non-deconvolved ima ge, and since they 
have been identified in previous studies (IBock et aL 2000; 
ITomono et al ] 1200 ill Table gives the positions for these 
sources and the 12.8 [im filter flux densities through aper- 
tures of 0.6 arcsec diameter (corresponding to the circles on 
Fig.EJ. 



4 DISCUSSION 

In the following, we compare the MIR maps to the adap- 
tive o ptics 3.5 urn map (Fig. 07) from NACO llRouan et alJ 
120041) . and to the HST/WFPC2 [OIII] A5007 map (Fig.08), 
obtained from the HST archive. The position ing of these 
maps is made according to lGalliano et al1ll2003l) : at 3.5 um , 
the co re coincides with the central engine (= radio source 
SI in iGallimore et allll996l) . and in [OIII], the NLR cloud 
B is located 0.1 arcsec north of the central engine. 

4.1 The core 

On the VISIR 12.8 ^.m image, we measure the following sizes 
for the core: on the non- deconvolved image, the N-S and E- 
W FWHM are respectively 0.44 arcsec and 0.29 arcsec, while 
on the deconvolved image, these values are 0.29 arcsec and 
0.1 9 arcsec. This is co nsistent with the FW HM derived 
bv | Tqmono et alJ (1200 ll) : 0.3 arcsecxO. 19 arcsec. iBock et alJ 
(2000) did not resolve t he core in the E-W direction, even 
at 0.1 arcsec resolution. iRouan et al.l (120041) give the N-S 
FWHM of 0.122 arcsec for the core at 3.5 urn. Their claim 
to have resolved the core in the E-W direction at 2.2 u.m has 
not been confirmed by the K-band visibi li ty me asurements 
with VLTI/ VINCI bv IWittkowski et alJ d2004l) . who find 
that the 2.2 u.m core comes from a region smaller than 5 mas 
at PA~45°. 

The comparison of the IR maps with the [OIII] map can 
help understand the N-S extension of the IR core. The very 
good spatial coincidence between the MIR sources and the 
NLR [OIII] emittin g clouds is well illust rated in Fig. As 
already discussed in lGalliano et all 1120031) . this suggests that 
the [OIII] emitting clouds are also MIR emitters. If it is so, 
the MIR emission in the central 0.4 arcsec diameter region 
is composed of the emission from the torus itself (the core), 
plus the emission of the NLR cloud B, 0.1 arcsec north, plus 
the (weaker) emission of NLR cloud A, 0.1 arcsec south. The 
MIR emission from the core appears then naturally extended 
along the N-S direction. This implies that the N-S extension 



is not likely to be related to any intrinsic property of the 
dusty torus. 

4.2 The extended emission 

The MIR extended emission displays an overall elongation 
along the NE-SW direction, similar to that of the NLR. Yet, 
a close look at the image reveals that the PA of the isophotes 
(drawn as ticks on Fig. 01) varies with their radius: the in- 
nermost isophote (r ~ 20 pc) is at PA=-4°, while the outer- 
most one (r ~ 300 pc) is at PA=+31°. The twisting of the 
MIR emission pattern, from PA=-4° to 31° between radii 
of 20 pc to 300 pc, highlights the presence of an inner spiral 
pattern which was already recognised on the 3.5 u.m image 
(Fig. 07). Interestingly, close to the central engine, the PA 
of the radio emission isophote (PA=-7° b etween the radio 
sourc es SI and S2, S2 located to the south. IGallimore et alJ 
1996) is very close to the PA value found on the MIR image 
at a scale of 20 pc. 

The [OIII] A5007 map (Fig 08) can be compared to 
the continuum subtracted [Nell] line image (Fig 04). Since 
this image is very similar to the 12.8 \na deconvolved image, 
but has a slightly lower angular resolution, we prefer to use 
the 12.8 U-m filter image for the comparison. The coincidence 
between the NE MIR sources and the NLR clouds is very 
good: the MIR core includes both the NLR clouds A and B, 
source NE1 corresponds to the NLR clouds C and E, NE2 
relates to the NLR cloud F, NE3 to the NLR cloud D and 
finall y NE4 to the N LR cloud G. As mentioned already in 
Galli ano et alJ 1120031) . this strongly supports the idea that 
all the MIR emission which is not related with the torus 
arises from dust associated with the NLR clouds. To the 
south, only weak [OIII] emission is found in coincidence with 
SW1 and SW2. 

The fact that in the [OIII] line the southern sources 
are weak compared to the northern ones, whereas in the 
[Nell] line they appear with comparable strengths, mainly 
occurs for an instrumental reason: the F502N filter ex- 
tends, in the rest frame of NGC1068, from -1400 kms -1 up 
to +200 kms"'. The HST/STIS spectra presented in 
ICecil et alJ J2002I) indicate that the northern clouds are 
moving towards us at ~ 500 to 800 km s _1 (remaining within 
the filter), while the southern [OIII] emission clouds are re- 
ceding with velocities around 500 kms -1 (falling outside the 
filter) . This implies that the southern part of the NLR must 
actually be brighter than it appears on the HST [OIII] map. 
Hence, the southern MIR sources SW1 and SW2 are also 
likely associated with NLR region clouds. In that respect, 
co mparison of the V I SIR 1 2.8 urn map with images obtained 
bv lThompson et alJ ll200 J) in the Ha and Pao? lines is en- 
lightening. There is an overall excellent agreement between 
the 12.8 um map and the Paa image: on the latter, in the 
SW quadrant two Paa emission knots can be seen, at dis- 
tances of about 2.2 arcsec and 3.4 arcsec from the central 
core. They match perfectly the source SW3 and the com- 
plex SW4 plus SW5. On the contrary, these knots are barely 
visible on the Ha image. This can be interpreted as extinc- 
tion of the southern cone of the NLR by a large-scale disc, 
roughly perpendicular to the NRL ax is. The extinctions de- 
rived bv lKraemer fc Crenshaw! (I2000T) are E(B - V) = 0.35 
for the southern cone of the NLR, and E(B - V) = 0.22 for 
the northern cone, supporting this hypothesis. 
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Using the 18.5 urn map jTomono et al.ll200ll. Fig. 09), 
we can derive the 12.8 um/18.5 um colour temperature for 
each aperture drawn in Fig. Q Assuming dust black-body 
radiation, we find that the temperature varies from 350 K 
to 200 K as the distance from the clouds to the central en- 
gine increases, up to ~4 arcsec. Hence, as this increases, 
the clouds become weaker in the MIR, and have cooler 
18.5 um/12.8 um colour, in agreement with the idea that 
the clouds are heated by the central engine. 

4.3 Could the MIR clouds be associated with 
young massive star clusters? 

The NLR clouds are generally thought to be photoionised 
by hard UV photons from the central engine. As we have 
just seen, the good correlation observed between the NIR- 
MIR sources and the NLR clouds supports the idea that 
the dominant contribution to the NIR-MIR emission comes 
from heated dust, surviving in the backside of dense clouds. 
Moreover, the spatial evolution of the MIR intensity and 
colours suggest heating by the central engine. Still, the ques- 
tion whether these clouds shelter star clusters has not been 
addressed yet. 

The observed size of the MIR sources, of the order of 
14 pc or less, is reminiscent of the size of Young Massive 
Star Clusters (YMCs) which are often found in starburst 
regions. The measured N-band fluxes and upper limits of 
the L-band fluxes of the identified NLR clouds (several Jy 
in the N-band, and several tens of mjy in the L-band) are at 
least one order of ma gnitude larger than t he typical fluxes 
of embedded YMCs llGalliano etai] 12001 and references 
therein). This means that the observed MIR flux comes from 
dust mainly heated by the AGN central engine, but does not 
discard the possibility for the NLR clouds to shelter YMCs. 

If the MIR sources were indeed hosting YMCs, one 
would expect important radio emission from the HII regions 
surrounding their massive stars or from their supernovae 
remnants. The lack of small-scale correlation between the 
radio and MIR emission argues against such a scenario, un- 
less the putative YMCs have a peculiar IMF with an upper 
mass cutoff lower than usually considered or unless the clus- 
ters are old enough for the massive stars to have already 
disappeared, hence suppressing the HII region emission and 
the SNR radio emission. 

Let us notice finally that the observed inner spiral pat- 
tern prompts also the question of the presence of YMCs in 
the NLR, as starburst activity is often induced by bars or 
mini-bars in the central regions of galaxies. So far, the argu- 
ments do not allow to draw any clear conclusion. In order to 
investigate whether YMCs are sheltered by the NLR clouds, 
the next step is to obtain NIR-MIR spectra at high spatial 
and spectral resolutions of the various clouds identified in 
this study. 
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